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ABSTRACT 

Context. The hot gas that fills the space between galaxies in clusters is rich in metals. In their large potential wells, galaxy clusters 
accumulate metals over the whole cluster history and hence they retain important information on cluster formation and evolution. 
Aims. We use a sample of 5 cool core clusters to study the distribution of metals in the ICM. We investigate whether the X-ray 
observations yield good estimates for the metal mass and whether the heavy elements abundances are consistent with a certain 
relative fraction of SN la to SNCC. 

Methods. We derive detailed metallicity maps of the clusters from XMM - Newton observations and we use them as a measure for 
the metal mass in the ICM. We determine radial profiles for several elements and using population synthesis and chemical enrichment 
models, we study the agreement between the measured abundances and the theoretical yields. 

Results. We show that even in relaxed clusters the distribution of metals show a lot of inhomogeneities. Using metal maps usually 
gives a metal mass 10-30% higher than the metal mass computed using a single extraction region, hence it is expected that most 
previous metal mass determination have underestimated metal mass. The abundance ratio of of-elements to Fe, even in the central 
parts of clusters, are consistent with an enrichment due to the combination of SN la and SNCC. 

Key words. X-rays: galaxies: clusters galaxies: clusters: general supernovae: general galaxies: abundances 



1. Introduction 

Since the first X-ray observations of the 7 keV iron line fea- 
ture in the 1970's bv iMitchell et alJ (Il976h we know that the 
[ intra-cluster medium (ICM) does not only contain primordial 
. elements but also heavy elements. As heavy elements are only 
produced in stars which reside mainly in galaxies the enriched 
material must have been ejected into the ICM by the member 
galaxies. Due to the large potential wells of galaxy clusters they 
retain all the enriched material, so it makes them excellent lab- 
I oratories for the study of nucleosynthesis and of the chemical 
enrichment history of the universe. Because the gas transfer af- 
fects the galaxy and galaxy cluster evolution, it is important to 
know when and how the enrichment takes place. 
' The components in a galaxy cluster interact with each other in 
' many different ways, thus to study the distribution of the ejected 
metals can give us important information on the mechanisms that 
transported the enriched gas into the ICM. 
Several processes were proposed to explai n the observed enrich - 
ment in the ICM: ram-pressure stripping ('Gunn & Gottlll972b . 
galactic winds (De Young 1978), galaxy-galaxy interactions 
(lGnedinlll99 8). AGN outflows (De Youns 1986; Hasinger_etJ] 
I2002h . intra-cluster supernovae (G erhard et al. 2002) and oth- 
ers. Simulations show an inhomogeneous dist ribution of the 
metals independent on the enrichment processes (iKapferer et al.l 
[200 6. Schindler & Diaferio 2008). Although AGN outflows as 
well as galaxy-galaxy interac tions can a dd metals to the ICM 
dKapferer et al.. .2005i: iHeath eTaP l20Q7b . simulations suggest 
that the metal enrichment of the ICM is primarily due to galac- 
tic winds and ram-pressure stripping. A detailed comparison be- 
tween the enrichment due to galactic winds and ram-pressure 
stripping revealed that these two processes yield different metal 
distributions and a different time dependence of the enrichment 



("Kapferer^etaD 12007b'). In massive clusters ram-pressure strip- 
ping provides a much more centrally concentrated distribution 
than galactic winds, because galactic winds can be suppressed 
in the cluster center while ram-pressure stripping is most effi- 
cient there due to the fact that the ICM density as well as the 
galaxies velocities are larger in the cluster center (IKapferer et al.l 
2006). 

X-ray spectra are the only measure for the metallicity of the 
ICM. The metallicity is derived mainly by measuring the equiv- 
alent width of the iron line once the continuum (almost en- 
tirely given by thermal bremsstrahlung) is known. With the 
first generation of satellites it was just possible to determine 
the r adial metallicity profile s (e.g. iMushotzky & LoewensteinI 
Grandi etaP 12004 ^. With deep observations of bright 
clusters of galaxies by Chandra and XMM-Newton satel- 
lites it is now possible to extract metallicities in certain 
regions of a galaxy clus ter and cons t ruct X-ray weighted 



metallicity maps (Schmidt 
O' Sullivan 

Sanders & Fabian 



Durret et al. 2005 

Werner etal. 2006; Sanders & Fabian 2006: 
20061: ISimionescu et al.ll2009l:lLovisari et al. 2009). 



t et al.1 120021 : ISanders et aP l2004l : 
m et al. 2005; Sauvaseot et al.'2005'; 

Hayakawa et alJ 



In this paper we present the results of the analysis of a sam- 
ple of 5 cool core clusters (Centaurus, Hydra A, Sersic 159-03, 
A496 and A2029) observed with XMM-Newton . Our first goal 
is to show that even in relaxed clusters the distribution of metals 
shows a lot of inhomogeneities that cause an underestimation of 
the metal mass. 

The paper is structured as follows: in Sect. 2 we present the data 
sets and data reduction techniques; in Sect. 3 we present spa- 
tially resolved measurements of metals abundances; in Sect. 4 
we present the metallicity and temperature maps; in Sect. 5 we 
determine the total, gas and metal mass of the clusters. A sum- 
mary of our conclusion is given in Sect 6. Throughout the pa- 
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per we assume Ho=70 km s"^ Mpc"\ Qa=0.73 and Dm=0.27. 
The elemental abundances presente d in this paper are given rel- 
ative to the solar abundances from I Anders & Grevesse (1989). 
We plot and tabulate values with error quoted at the 90% level 
of confidence. 

2. X-ray observations and analysis 

2. 1. Sample selection 

Our first aim is to derive the metallicity maps of the cluster. To 
determine metallicities in many different regions of a cluster a 
lot of photons are required. Therefore we selected clusters ob- 
served with sufficiently long exposure from the XMM -Newton 
archive. We found 5 clusters suitable for our study: Centaurus, 
A496, Sersic 159-03, Hydra A and A2029. We did not take the 
Perseus Cluster into account, because in literature there is al- 
ready a detaile d metallicity map obtained using Chandra data 
(ISanders et al.ll2b 05). We decided to include Centaurus, because 
ISanders & FabianI (|2006) derived abundance maps for several el- 
ements but only for the inner part of the cluster (-80 kpc.) 

2.2. Data reduction 

Observation data files (ODFs) were retrieved from the XMM 
archive and reprocessed with the XMM-Newton Science 
Analysis System (SAS) vT.l.O. We used tasks emchain and 
epchain to generate calibrated event files from raw data. 
Throughout this analysis single pixel events for the pn data 
(PATTERN 0) are selected, while for the MOS data sets the 
PATTERNS 0-12 are used. In addition, for all cameras events 
next to CCD edges and next to bad pixels were excluded 
(FLAG==0). 

The data were cleaned for periods of high background due to 
the soft proton solar flares using a two stage filtering process. 
We first accumulated in 100 s bins the light curve in the [10-12] 
keV band for MOS and [12-14] keV for pn, where the emis- 
sion is dominated by the particle-induced background, and ex- 
clude all the intervals of exposure time having a count rate that 
dev iated by more than Scr from the mean (see iPratt & ArnaudI 
12002) for details. After filtering using the good time intervals 
from this screening, the event lists was then re-filtered in a sec- 
ond pass as a sa fety check for possible flares with soft spectra 
dN^alainen et al. 2005; Pradas & Kerp 2005). In this case fight 
curves were made with 10 s bins in the full [0.3-10] keV band. 
The resulting exposure times after cleaning are listed in Table [T] 
For the background subtraction we used a combination of 
blank-sky maps and c losed-filter observations as done by 
ISimionescu et al.l (l2009l) . From deep sky observations collected 
with XMM, we selected the data with the most similar back- 
ground for each cluster. Both the blank-sky and the closed fil- 
ter events were selected by applying the same PATTERN se- 
lection, vignetting correction, flare rejection criteria and point 
source removal used for the observation events. In addition we 
transformed the coordinates of the background files such that 
they were the same as for the associated cluster data set. We cal- 
culated the count rates in the hard energy band (10-12 keV for 
MOS and 12-14 keV for pn) outside of field of view (OOFOV) 
for each observation, blank sky maps and closed filter observa- 
tions. For each detector we added to the corresponding blank sky 
background set a fraction of the closed filter observation in order 
to compensate for the difference between the OOFOV hard-band 
count rate in the observation and in the blank-sky data. 
To correct for the vignetting effect, we used the photon weight- 



Table 1. Cluster sample, sorted in order of increasing redshift. 
The exposure column indicates the net exposure time, after ex- 
cluding time intervals of strong background flares. 



Cluster 


redshift Orbit Rev. 




exposure (ks) 










MOSl 


MOSl 


pn 


Centaurus 


0.0114 


0379 


28.0 


35.8 


31.5 






1213 


107.5 


106.5 


80.6 


A496 


0.0329 


1405 


58.8 


58.9 


41.8 






1501 


61.2 


60.0 


46.4 


HydraA 


0.0539 


1359 


81.6 


82.4 


55.1 


Sersic 159-03 


0.0564 


0540 


92.9 


94.2 


73.9 


A2029 


0.0773 


0496 


12.1 


12.1 


9.2 






1576 


34.3 


37.6 


19.5 






1577 


41.7 


43.8 


27.3 






1578 


33.5 


26.4 


18.6 



ing method (lArnaud et al .112001"). The weight coefficients were 
computing by applying the SAS task evigweight to each event 
file. Point sources were detected using the task ewavelet in the 
energy band [0.3-10] keV and checked by eye on images gen- 
erated for each detector. We produced a list of selected point 
sources from all available detectors and the events in the corre- 
sponding regions were removed from both the blank field and 
the observation data set. 

2.3. Spectral analysis 

All the spectral fits were perf ormed with the XSPEC package 
(version 12.5.0, Arnaud 1996). In order to model the emission 
from a single (or multi) tempe rature plasrna we fit the spectra 
with an APEC (+APEC) model (ISmith et al.ll200lh multiplied by 
the Galactic hydrogen column density, N//, fixed at the Galactic 
values |Dickev & Lockman 1990) through the wabs absorption 
model (IMorrison & Mc Cammon 1983). To determine the abun- 
dance profiles a VAPEC(+VAPEC) is used. We fit jointly MOSl, 
M0S2 and pn spectra, enforcing the same normalization value 
for MOS spectra and allowing the pn spectrum to have a sepa- 
rate normalization. In the spectral fitting we used the 0.3-10 keV 
energy range. Because of the low number of counts and the high 
background at large radii the spectra are fitted in the 0.3-9 keV 
band and 0.5-7.5 keV energy range for MOS and pn respectively. 
For the pn data we excluded the energy above 7.5 keV in order 
to discard strong instrumental lines around 8 keV. The redistri- 
bution and ancillary files (RMF and ARE) were created with the 
SAS tasks rmfgen and arfgen for each camera and each region 
that we analyzed. 

3. Abundances Profiles 

In order to know how the ICM has been enriched, we need to 
measure the amount and distribution of metals in the ICM. Since 
Ar and Fe are synthetized mainly in SNe la, O in core-collapse 
supernovae (SNCC), S and Si in both SN la and SNCC, the mea- 
sure of these elements can give hints on the past enrichment pro- 
cess in the ICM by supernovae. 

For each cluster we extracted the spectra from several annular 
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Fig. 1. Radial profiles of the (a) iron, (b) oxygen, (c) silicon, (d) sulfur and (e) argon in solar units, derived from the EPIC spectra. 



regions. Metal abundance profiles are determined for O, Si, S, 
Arand[Fe;_We fitted the data with the procedure presented in 
iLovi sari et al. I([2b09) to avoid the degeneracy of the parameters: 
(1) we fitted the data with an absorbed APEC model in the 0.4-7 
keV band to obtain the temperature (metallicity and normaliza- 
tion are considered free parameters); (2) we fixed the tempera- 
ture and use a VAPEC model in the same energy band to deter- 
mine the iron abundance (O, Mg, Si, S, Ar, Ni were left free. The 
other elements were fixed to the solar value); (3) we kept tem- 
perature and iron fixed to measure the oxygen abundance in the 
0.4-1.5 keV band; (4) we fix the values of temperature, iron and 
oxygen to estimate the silicon, sulfur and argon abundances in 
the 1.5-5 keV band. Fixing to the solar value the abundances of 
elements, which cannot be significantly detected, may introduce 
a bias in the abundance measurements of other elements if the 
average cluster abundance is not solar. Thus, as a safety check 
we fitted the spectra with C, N and Ne fixed at 0.3. Apart from 
Centaurus the resulting abundance of the interested elements (O, 
Si, S, Fe) do not change significantly. Furthermore, we note that 
in general, when we examined all abundances to be free in the 
fits, the resulting parameters did not change within the statistical 
errors. 

Since there are still cross-calibration issues at energies below 
~ 2keV we fit the MOS and pn spectra separately to investigate 
how robust the derived values are. In general the values appear 
to be consistent within the errors in the two instruments although 
for some annuli the abundances of several elements are clearly 
not. In particular O shows a strong discrepancy for Centaurus, 
A2029 and Hydra A up to 0.5 arcmin while Si is not consistent 



for Sersic and Hydra A for the annulus between 0.5 and 1 ar- 
cmin. The systematic difference between MOS and pn is 3% for 
Fe, 17% for O and 11% for Si. 

In general the single temperature model (IT) provided a 
good fit to the spectra. However, in the inner region of 
all the clusters a 2T model significantly improved the fit 
with res pect to a sin gle temperature m odel as al ready found 
by Sanders & FabianI (12006) and Matsushita et al.l ([2007a) for 
Centaurus, b y Isimionescu et al. I (l2009ir for Hydra A and by 
Ide Plaa et all ('2006') for Sersic 159-03. The Fe abundances 
changed significantly as compared to the single temperature 
fit (the so-called Fe-bias, se e for example iBuotd 120001 and 
iMolendi & Gastaldellol [200lh . For these regions we use the 
abundances derived from this 2T model for the discussion be- 
low. 

The abundance value for Si, S, Ar and Fe and their radial varia- 
tion look quite similar to each other. They commonly decline to 
about 1/4 of the central value in the outermost annulus (see Fig. 
[T]). The error on the Si, S, and Ar abundance determination is 
larger than on iron. Nevertheless, their profiles indicate a central 
increase similar to the iron distribution. 

The O profile looks more complex compared with the other ele- 
ments. It shows several discontinuities, probably due to the fact 
that O is strongly related to episodes of star formation, in fact 
the O abundance seems to increase in the outer rings, where we 
expect higher star formation. On the other hand, while in the 
central regions the flux is high enough to get an accurate mea- 
surement, at large radii, the oxygen abundance could be overes- 
timated because the oxygen in the galactic foreground emission 
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Table 2. Fit results for spatially resolved EPIC spectra. For each cluster we show also the results obtained fitting the spectra of a 
radius encompassing an overdensity of ~ 4500. 



Target 


radius 


k1 1 


KlZ 


re 


(J/re 


5l/re 


5/re 


Ar/Fe 


X Id.o.f. 




arcmin 


keV 


keV 














Centaurus 


0-0.25 


1.763+0.015 


0.778+0.002 


1.190+0.002 


0.48+0.03 


0.97+0.03 


1.01+0.05 


0.77+0.12 


1492/984 




0.25-0.5 


1.853+0.008 


0.846+0.003 


1.276+0.017 


0.49+0.04 


1.03+0.03 


0.98+0.04 


1.05+0.10 


2149/1290 




0.5-0.7 


1.969+0.006 


1.001+0.010 


1.336+0.023 


0.58+0.04 


1.11+0.05 


1.07+0.06 


0.88+0.10 


1917/1303 




0.7-1 


2.120+0.010 


1.078+0.009 


1.360+0.017 


0.56+0.04 


1.14+0.04 


1.02+0.04 


0.80+0.09 


2390/1510 




1-1.4 


2.462+0.009 


1.275+0.010 


1.500+0.015 


0.44+0.03 


1.02+0.04 


1.02+0.04 


0.77+0.08 


2652/1743 




1.4-2 


2.755+0.008 




1.416+0.010 


0.39+0.02 


1.07+0.03 


0.98+0.04 


0.69+0.06 


3168/1963 




2-2.8 


2.998+0.008 




1.269+0.010 


0.42+0.03 


1.09+0.03 


0.95+0.03 


0.76+0.07 


3722/3110 




2.8-4 


3.368+0.008 




0.932+0.008 


0.50+0.04 


1.10+0.05 


0.92+0.05 


0.72+0.10 


3363/2258 




4-5.6 


3.685+0.012 




0.683+0.008 


0.56+0.05 


1.11+0.06 


0.79+0.07 


0.89+0.15 


3222/2448 




5.6-8 


3.745+0.013 




0.504+0.008 


0.81+0.08 


1.15+0.08 


0.85+0.10 


0.79+0.22 


3159/2386 




8-13 


3.800+0.016 




0.434+0.010 


0.92+0.09 


1.27+0.09 


0.94+0.09 


0.76+0.43 


3503/2436 




0-13 


3.773+0.016 


1.662+0.007 


0.649+0.010 


0.64+0.04 


1.30+0.07 


1.13+0.08 


1.01+0.15 


2590/2051 


A496 


0-0.5 


2.903+0.022 


1.332+0.012 


0.757+0.012 


0.61+0.05 


1.25+0.06 


1.23+0.08 


0.90+0.17 


2115/1688 




0.5-1 


3.675+0.071 


1.710+0.027 


0.626+0.009 


0.71+0.05 


1.27+0.06 


1.24+0.07 


1.07+0.18 


2347/1951 




1-2 


4.462+0.140 


2.118+0.018 


0.524+0.007 


0.76+0.07 


1.32+0.07 


0.95+0.09 


0.95+0.21 


2897/2263 




2-3 


3.944+0.018 




0.464+0.009 


0.93+0.13 


1.37+0.12 


0.99+0.13 


0.91+0.28 


2734/2204 




3-4 


4.060+0.022 




0.398+0.009 


1.08+0.15 


1.40+0.13 


0.83+0.17 


1.05+0.39 


2491/2114 




4-6 


4.100+0.025 




0.340+0.009 


1.20+0.18 


1.44+0.20 


0.76+0.15 


0.91+0.39 


2674/2095 




6-9 


3.981+0.025 




0.246+0.010 


1.54+0.28 


1.59+0.23 


0.69+0.24 


1.10+0.34 


2700/2026 




0-5 


3.453+0.024 


2.081+0.018 


0.525+0.010 


0.71+0.11 


1.16+0.12 


0.83+0.13 


0.88+0.24 


2760/2214 


Sersic 159-03 


0-0.5 


2.783+0.027 


1.937+0.013 


0.460+0.012 


0.67+0.06 


1.08+0.07 


1.06+0.12 


0.74+0.24 


1631/1429 




0.5-1 


3.103+0.022 


2.192+0.026 


0.433+0.008 


0.67+0.08 


1.15+0.10 


1.15+0.10 


0.83+0.23 


1731/1432 




1-2 


2.574+0.012 




0.330+0.006 


0.74+0.09 


1.42+0.11 


0.59+0.12 


0.83+0.28 


1846/1587 




2-3 


2.505+0.022 




0.250+0.010 


0.88+0.20 


1.64+0.23 


1.40+0.27 


1.16+0.63 


1386/1221 




3-4 


2.357+0.032 




0.178+0.011 


0.94+0.18 


1.97+0.19 


1.35+0.32 


0.77+0.54 


1226/1007 




4-6 


2.063+0.020 




0.133+0.010 


3.00+0.20 


2.33+0.29 


2.07+0.47 


1.90+1.60 


1201/1012 




0-2.5 


2.590+0.018 


1.303+0.026 


0.405+0.007 


0.76+0.06 


1.09+0.05 


0.83+0.05 


0.79+0.14 


2587/2030 


Hydra A 


0-0.5 


5.064+0.015 


2.232+0.077 


0.488+0.009 


0.94+0.11 


1.23+0.09 


1.00+0.13 


1.25+0.30 


1979/1730 




0.5-1 


5.354+0.015 


2.423+0.245 


0.413+0.010 


0.85+0.12 


1.40+0.10 


1.11+0.15 


1.15+0.36 


2029/1796 




1-2 


5.027+0.022 


2.165+0.230 


0.355+0.009 


0.88+0.13 


1.33+0.15 


1.00+0.17 


0.39+0.81 


2271/1906 




2-3 


5.077+0.020 


2.267+0.153 


0.269+0.010 


0.74+0.19 


1.37+0.20 


1.31+0.24 


1.08+0.47 


1929/1761 




3-4 


3.879+0.035 




0.273+0.013 


0.84+0.28 


1.13+0.28 


0.94+0.31 


0.73+0.56 


1790/1657 




4-6 


3.760+0.034 




0.329+0.019 


0.88+0.31 


1.31+0.40 


1.06+0.36 


0.00+0.16 


2076/1789 




6-8 


3.706+0.090 




0.333+0.020 


0.31+0.37 


0.97+0.45 


0.27+0.27 


0.99+0.77 


2132/1724 




0-3 


5.090+0.127 


2.183+0.039 


0.370+0.008 


0.66+0.06 


1.21+0.07 


1.02+0.07 


0.66+0.16 


3006/2505 


A2029 


0-0.25 


7.715+0.040 


2.891+0.127 


0.744+0.013 


0.93+0.11 


1.22+0.18 


0.95+0.20 


0.06+0.33 


1750/1657 




0.25-0.5 


8.415+0.041 


3.394+0.142 


0.689+0.013 


0.90+0.12 


1.15+0.16 


1.01+0.21 


0.00+0.17 


2461/2112 




0.5-0.7 


8.190+0.047 


3.214+0.167 


0.568+0.016 


0.92+0.17 


0.98+0.20 


0.99+0.29 


1.46+0.65 


2288/2054 




0.7-1 


7.886+0.047 


2,751+0.182 


0.516+0.012 


0.92+0.18 


0.76+0.21 


0.97+0.30 


0.31+0.55 


2467/2220 




1-1.4 


8.191+0.135 


2.068+0.140 


0.402+0.011 


0.78+0.24 


0.74+0.28 


0.95+0.36 


0.00+0.26 


2508/2303 




1.4-2 


8.415+0.065 


2.435+0.187 


0.400+0.012 


1.02+0.24 


1.05+0.29 


0.65+0.36 


1.34+0.82 


2571/2377 




2-2.8 


7.388+0.060 




0.366+0.015 


2.21+0.31 


0.55+0.34 


0.88+0.34 


0.00+0.41 


2673/2303 




2.8-4 


7.704+0.066 




0.331+0.011 


0.79+0.37 


0.76+0.40 


0.48+0.36 


0.00+0.16 


2682/2299 




4-6 


8.215+0.085 




0.341+0.016 


2.43+0.49 


0.15+0.34 


0.97+0.60 


0.00+0.58 


3286/2470 




0-3 


8.100+0.087 


2.385+0.123 


0.470+0.016 


1.12+0.13 


0.94+0.13 


1.00+0.16 


0.43+0.07 


3490/2791 
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starts to play an important role. This effect is most apparent if the 
background is high with respect to the source emission. Apart 
from Centaurus, in the outermost bin of each cluster we found 
approximately a S/N of 4, thus the O abundance should be used 
with caution there. 

We computed abundance ratio of O, Si, S and Ar over Fe as a 
function of the projected radius (see Table O. The ratios Si/Fe 
and S/Fe are consistent with a constant values around 1-2 and 
0.7-1.5 respectively, while the 0/Fe ratio for the innermost re- 
gion is lower around 0.5-1 (see Fig. O. In addition the 0/Fe 
ratio suggests some increase with radius with the exception of 
Hydra A for which both 0/Fe and Si/Fe ratio seems to be con- 
stant. If we do not consider the outermost bin the 0/Fe show a 
sli ght increase with radius i n agreement with the results obtained 
bv lSimionescu et al.1 (l2009l) . 

Using Suzaku o bser vations. dSato et al.1 l2007l 12008 allbllch . 
iTokoi et al.1 ('2008) and [Matsushita et al.1 (l2007bb have presented 
abundances of groups and clusters of galaxies. All systems 
show very similar value o f Si/Fe ratio, to be 1-1.5, in good 
agreement with our results. iTamura et al.l (|2004|) reported abun- 
dance ratios for 19 clusters (among them A496, Hydra A, 
Sersic but with shorter exposure time) studied with XMM- 
Newton , and the mean Si/Fe was -1.4. Their 0/Fe, -0.6, is 
bit lower than our results. A ratio of 1-2 for the Si/Fe ratio was 
fou nd also b y de Grandi & Molencii| (l2009h analyzing 26 clus- 
ters, ^^^^s & Fabian (2006) showed the abundances ratios for 
the Centaurus cluster with Chandra and XMM-Newton , and the 
radial abundance ratios of 0/Fe and Si/Fe to be 0.5-1 and 1-1.5 
respectively, were consistent with our results. 
Therefore, in general our results are in agreement with previous 
studies suggesting that cluster and groups have passed the same 
metal enrichment process in the ICM. 

3. 1. Number ratio of SNCC to SN la 

We investigated the SNe la and SNCC (including type lb and 
Ic) contr ibution to the I CM m etals following the method pre- 
sented in'LovisarietaLl (l2009h . Using the SNe la of the WDD2 
model (t he curren tly favou red SN la explosion s cenario) adapted 
from Iw amoto et all (1 1999b and SNCC yields bv lTsuiimoto et al.l 
(1995) we computed which SNCC/Ia ratio reproduces the ob- 
served abundances better. We note that although the models 
adopted here (SNe yields, Salpeter IMF, spherical symmetry in 
the ICM, etc.) are probably too simplified and we should take 
into account other effects (i.e metals locked in stars), we want to 
show if there is agreement between the measured abundances 
and the theoretical yields for different elements. Table [3] and 
Fig. [3] summarize the results. We found good agreement in the 
O, Si and Fe abundances for A496, Sersic, and Centaurus for 
which the relative fraction of SNCC seems to increase toward 
the outskirts, and for Hydra A for which the relative fraction of 
SNCC seems to be constant. Due to the large error bars in the 
SN determination of A2029 it is difficult to say whether there 
is an agreement between the same elements. The flatness pro- 
file of the relative fraction of SNCC in Hydra A up to 370 kpc 
could be a c onsequence of the mixing due to the central AGN. 
ISimionescu ^et al. ( 2009) showed that the 0/Fe in the cooler gas 
extending in arm-like structures caused by the central AGN is 
consistent with the average 0/Fe ratio in the inner 3\ Concerning 
A2029, we found an agreement between the O, Si, and Fe only 
in a few regions probably due to the fact that the spectral lines 
are weaker at high temperature and the determination of the el- 
ements becomes difficult. We note that when the S/Fe ratio in- 
creases the SNCC fraction computed using Si/Fe ratio is not in 



Table 3. Number ratio of SNCC to SN la (N///N and integrated 
number of SNCC computed using the abundance of O for each 
considered region. For each cluster we show also the results ob- 
taining fitting the spectra of a radius encompassing an overden- 
sity of -4500. 



Cluster 



region 



0/Fe 



N///N;, 

Si/Fe S/Fe 



Centaurus 



A496 



Sersic 



Hydra A 



A2029 



0-0.25 
0.25-0.5 
0.5-0.7 
0.7-1 

1- 1.4 
1.4-2 

2- 2.8 
2.8-4 
4-5.6 
5.6-8 
8-13 
0-13 



0-0.5 
0.5-1 

1- 2 

2- 3 

3- 4 

4- 6 
6-9 
0-5 



0-0.5 
0-5-1 

1- 2 

2- 3 

3- 4 

4- 6 
0-2.5 



0-0.5 
0.5-1 

1- 2 

2- 3 

3- 4 

4- 6 
6-8 
0-3 



0-0.25 
0.25-0.5 
0.5-0.7 
0.7-1 

1- 1.4 
1.4-2 

2- 2.8 
2.8-4 

4-6 
0-3 



1 09+008 

1 ll+O-^l 

1 30+^-^^ 

93+007 

1 14+0.11 

^•^^-0.11 

1 30+^-^^ 
9 44+0.34 

1 C4+O.IO 



77+012 

^•^^-0.14 

1 07+0.24 
^••^ -0.23 
1 C9+O.I8 
^•^^-0.20 

97+016 

1 09+0.25 

-"^••^^-0.22 
1 07+0.19 

-0.35 
9 17+0.53 

^' ^ ' -0.48 
9 00+0.37 
^••^•^-0.37 



95+^-^2 

^•^-^-0.34 
^' ^^-0.25 
^•^^-0.46 

1 03+^-^^ 
1 03+^-^^ 
74+0-29 

^•^^-0.25 

34+0-34 



<0.53 
<1.26 

1 Q9+0.78 



1 oc+0.14 
^"^^-0.13 

1 75+0-21 

^' ^^-0.17 
1 oc+0.23 
i.8^_Q2i 

9 OO+0.48 

3 03+^-^^ 

^•^^-0.59 

4C+O.54 
-^•^^-0.51 

^•^^-2.11 

1 79+0.36 

-0 . 31 



1 69+^-^^ 

^•"^-0.32 
9 04+0.52 
^•^^-0.43 
9 44+0.42 
^•^^-0.40 

2 70+0-76 

^' ^^-0.88 
9 Q4+0.99 
^•^^-0.81 

''•^ -0.91 

4 35+2-31 

^•^^-1.85 

1 61+^-^^ 
^•'^^_n.55 



1 07+0.55 

^••^^-0.56 

19+1.22 
•^•^ -0.90 

<0.54 

1 25+1-^^ 

^•^^-0.85 

<1.27 
<2.69 
<1.10 
<0.60 



1 61+0-^^ 

1 61+^-24 

9 17+0.26 
^' ^ ' -0.32 
9 97+0.76 
-0.63 

2 50+°-^^ 

^"^^-0.61 

25 0+^-^^ 
1 89+^-^^ 

^•^^-0.20 



1 14+0.31 
1 56+^-^^ 

^"^"-0.46 
+0.63 
0.58 



2.70 



4 54+2-12 

^•^^-2.42 

14 3+10-7 
^^••^-5.10 
1 90+0.23 

^•^^-0.23 



1 33+1-11 

^•^^0.0.87 
9 10 + 1.00 
^•^'^-0.86 

<0.88 

5-26!^:^^ 

4 Q<+6.76 



<0.19 



9 00+0.32 

^••^"-0.38 

9 97+0.50 
^•^^-0.42 
1 QO+0.62 
^•^^-0.55 
9 1 9+0.65 
^•^^-0.86 
9 00 +0.90 
^••^•^-0.81 

<1.61 

1 CQ+0.16 
^•^^-0.18 



1 61+^-^^ 

^•'^^-0.48 

2 70+0-63 

^•^^-0.66 

2 50+^-^^ 

^•^^-0.80 

2 70+1-^0 

^•^^-1.06 

1 41+1-^1 

'5Q + 1.62 
Z. Jo_j^ 22 

<3.12 

1 o<+0.37 
^•^^-034 



<0.86 

1 4Q+1.28 
^•^^-1.00 

<1.41 

3 84+"^-^^ 

•^•°^-2.44 

<2.63 



1 03+0-56 

^•^-^-0.58 



9 44+0.42 
^•^^-0.36 

QQ+0.40 

9 44+0.59 
^•^^-0.59 
9 44+0.68 
^•^^-0.60 

1 96+0-^2 

^•^'^-0.71 
9 70+0.93 
^•/^-0.82 
Q QQ+4.28 
^•^^-2.30 

1 QQ + 1.16 

12 5+12.5 

^^•^-4.80 
o 19+0.58 
•^•-^ -0.58 



1 89+1-^5 

^•^^-0.82 
1 C4+0.84 
-0.69 

0.81!^:?: 

<0.76 
<0.97 
<2.55 
<0.46 
<1.61 



66+^ 



<2.13 
<2.85 
<3.57 
<3.27 
<3.84 
<0.95 
<2.71 

<8.32 
<2.22 



Ar/Fe 



<1.67 

6 25+^-^"^ 

"•^-^-2.40 

2 50+2-04 

^•^^-1.54 

1 35+3-19 

^"^^-0.39 
^•^'J_0.84 

<0.63 

0-86!^:?? 

<1.61 

2 63+3-25 
<5.00 
<12.5 

5.oo;5-QQ 



<2.49 

12.5:24-7 

3.03!5-29 

3 03+3-22 

9 Q4+3.31 
-2.15 

4 50+12.1 

^•^^-3.80 

7 70+42 

^•^^-6.30 

<8.33 



<2.75 
<4.91 
<5.03 

<24.0 

<3.33 



19 c+10.0 
^^•^-10.0 



<19 

<24.5 

<4.76 
<1.61 



<2.50 
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Fig. 2. Radial profiles of the (a) 0/Fe , (b) Si/Fe and (c) S/Fe ratio in solar units. 



Centaurus 



A496 



A2029 




Hydra A 



Radius (arcmin) 



Fig. 3. Radial profiles of the relative fraction of SNCC using the 0/Fe (red) and Si/Fe (black) ratio. 



agreement with what obtained using 0/Fe ratio. Furthermore, 
even for most of the considered regions in the other clusters, 
the S/Fe ratio suggests a very low contribution to the enrichment 
due to SNCC that is in contradiction with the idea of an early en- 
richment due to massive stars. In other words, either the yields 
of sulfur from theoretical works are underestimated or when the 
abundances of Si is high t he determination of sulfu r becomes 
very difficult. We note that tea umgartner et al.ll2QQ5l) found that 
clusters with a temperature between 6 and 8 keV (as for A2029) 
show a general increasing of sulfur corresponding to a decreas- 
ing of silicon. 



Apart from A2029, the abundances of the other four clusters 
are consistent with a SNCC to SN la ratio of around 1.5-3. 
This ratio means that almost 75-85% of Fe, 15-25% of Si and 
< 20% of S were synthetized by SN la. Using deep qbserva- 
tio ns of 2A 0335 +096 and Sersic 159-03 "Werner et al." (20061) 
an d|de Plaa et al. 1|2006) found a ratio of 2-2.5 in agreement with 
our result. B a sed o n XMM-Newton observations of 22 clusters 
Ide Plaa et al.l (12007) found a ratio in the range 1.7-3.5 depend- 
ing on the supernovae models (using a WDD2 model, as in this 
work, they obtained 1.7 in good agrement with our result). 
The efficiency of the metal enrichment may depend on param- 
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eters as age, star-formation efficiency and contribution from cD 
galaxies. However, the relative contribution of SNe la and SNCC 
and the processes of metal mixing in the ICM seem t o be quite 
similar among different clusters as claimed also by iSato et alJ 
(l2008ah . Our results are therefore consistent with an enrich- 
ment due to the combination of SN la and SNCC. This con- 
clusion contrasts with the general views that the central iron 
excess in cool core clu sters is from Type la supernovae alone 
([Matsushita et al.ll2Q03l) . 

4. 2D Distribution of Metals 

In order to study the distribution of metals, we prepared adap- 
tively binned abundance maps of the clusters. Appropriately 
cleaned data sets for all the three cameras, with point sources re- 
moved were used to create the source spectra. To obtain a metal- 
licity measurement with a good accuracy requires a high statis- 
tic. Thus, to ensure an acceptable error also in the outskirts of the 
clusters we set a minimum count number (-5000 source counts 
per region) necessary for proceeding with the spectral fit. The 
spectral regi ons for the map were selected following the method 
presented in lLovisari et al.l (l2009l) that we can summarize in this 
way: a square region centered on the X-ray peak was defined to 
include the area with high surface brightness. The region size of 
the pixels was optimized to be as small as possible by splitting it 
into horizontal or vertical segments through its center, while in- 
cluding at least 5000 source counts. For all the selected regions, 
spectra were extracted for source and background in all three 
cameras. 

The obtained metallicity maps are shown in Fig. |4] and [5l The 
metallicity distribution appears very inhomogeneous for all the 
clusters. For Centaurus and Sersic 159-03, there is a peak in the 
center and then it decrease in the outskirts while A496, A2029 
and Hydra A show high metallicity clumps both in the center 
and in the outskirts. We note that, since Centaurus is at very low 
redshift, in these observations we are looking in the very cen- 
tral part of the cluster (r<200 kpc) compared with the others 4 
clusters, and it could explain its different shape. On the other 
hand, Sersic, for which we map the metal distribution for more 
than 350 kpc we observe the same shape of Centaurus. Since 
it is difficult to distinguish real metallicity clumps from statisti- 
cal noise we quantified the inhomogeneities through the signifi- 
cance maps. First, we smoothed the metallicity profiles applying 
a Savitzky-Golay filter and we subtracted it from the metallicity 
maps. Finally we divided each bin of the resultant maps by the 
uncertainty in the pixel metallicity. The results are shown in Fig 
[6l The blue and red spots represent regions that deviate significa- 
tively (99% c.l.) from the average profile. 
There are several maxima visible in the metal distribution, which 
are not associated with th e cluster center. From simulations 
(iKapferer et al.ll200"8[ 1200 9') we know that the maxima are typ- 
ically at places where galaxies just have lost a lot of gas due to 
the ram-pressure. Since the gas lost by galaxies is obviously not 
mixed immediately with the ICM at the place where we observe 
a metal blob we should observe also a low temperature due to 
the fact that the gas in galaxies is cooler than the ICM. Thus, 
we produced the temperature maps (fitting the spectra with a 
single temperature model) of the clusters with the same spatial 
resolution obtained for the metal maps and then we plotted the 
abundance of bins against their temperature, which are shown in 
Fig. [71 Since we are searching for cool high-metallicity clumps 
due to the ejection of gas from galaxies and not the cool high- 
metallicity bins found in the cool cores, we did not plot the inner 
bins (where we used a two temperature model to fit the spectra). 
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Apart from Centaurus, we see a deviation from the expected 
temperature-metallicity relation. The deviation could be due to 
the combination of the iron bias and inverse iron bias effects. 
The iron bias effects cause an underestimation of the metallicity 
when we pretend to fit with a single temperature model a plasma 
that is, ins tead characterized by a co mbination of different tem- 
peratures (lBuotdl200Q[ iRasia et al.l 12008). On the contrary the 
inverse iro n bias effect cause an o verestimation of the metal 
abundance (iGastaldello et al.ll2010l) . The combination of these 
two effects could explain the spread in the distribution of Fig. [71 
On the other hand these effects seems to be more important in 
the range of temperature between 2 and 4 keV, while our clus- 
ters show also regions with temperature higher than 4 keV. One 
possible explanation could be that the ej ected gas, with T < ] 
keV and m etallicity in the range 0.5-1.5 ( Matsushita et al.|[2000l 
lAthey & Bregman 2009), will be heated up to the temperature 
of the surrounding gas (ICM) on the shorter time- scale than that 
of metal mixing. In this case, after a while we should observe 
a region of high metallicity (not yet dispersed) and high tem- 
perature (heated up at the ICM temperature). Another possible 
explanation could be related to the number of intracluster super- 
novae. In fact during ram-pressure stripping events a lot of stars 
are forming in the tail of stripped gas. The stars are evolving 
and exploding as SNe directly in the ICM and they can enrich 
the ICM very efficiently. In this case we should see a clumps of 
high metallicity (due to SNe explosions) and high temperature. 
Obviously, more complex heating and cooling processes are at 
work, thus the simple picture of stripped gas does not hold. 

5. Mass determination 

Using simulated galaxy clusters iKapferer et al.l (l2007ah showed 
that the more inhomogeneous the metals are distributed within 
the cluster, the more underestimated is metal mass. They showed 
that the true metal mass in the inner parts (r<500 kpc) of galaxy 
cluster can be up to three times higher than the metal mass ob- 
tained by X-ray observations. They suggested that the discrep- 
ancies are due to the fact that the metallicity is not constant 
throughout the extraction area, thus the integration of thermal 
bremsstrahlung and of line emission can lead to underestimated 
metal masses. 

We used the metallicity maps to estimate the metal mass in the 
center parts of the clusters and compare the results with the esti- 
mations obtained from a single extraction area of the cluster. 
First, we computed a background subtracted, vignetting cor- 
rected, radial surface brightness profile in the 0.3-10 keV en- 
ergy band for each cluster. All of the X-ray point sources were 
escluded from the data. The annuli were chosen such that all 
the widths are larger than the FWHM of the point- spread func- 
tion (PSF) at that radius. With this choice, all the bins con- 
tain at least 2000 counts after background subtraction. With 
this very good statistics, the error bars of the surface bright- 
ness are very small. The profiles were fitted using a /5-model 
(Cavaliere & Fusco-Femianal976.) : 

/ ^2 \ -3/5+0.5 

S(r) = So\l^-^J (1) 

where is the core radius. The advantage of using a yS-model 
to parametrize surface brightness is that assuming hydrostatic 
equilibrium and spherical symmetry the gas density and total 
mass profile can be recovered analytically and expressed by the 
simple formula: 

ngas(r) = nod + x")-'^'^ (2) 
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Fig. 4. left: Metallicity maps based on spectra from all three EPIC camera, right: Temperature maps, obtainded fitting the spectra 
with a single temperature model, with the same resolution of the metal maps. The scale for the metallicity and temperature is in 
solar units and keV, respectively. 



L. Lovisari et aL: Inhomogeneous Metal Distribution in the Intra-Cluster Medium 



9 



-42 36 00 - 




23 14 40 23 14 20 23 14 00 23 13 40 23 13 20 

RIGHT ASCENSION (J2000) 



0.1 0.2 0.3 0.4 



H- -12 05 00 




09 18 40 


1 1 

09 18 20 09 18 00 

RIGHT ASCENSION fJ2000^ 


1 

09 1 7 40 


1 -^^^m 



23 14 40 23 14 20 23 14 00 23 13 40 

RIGHT ASCENSION (J2000) 




09 18 20 09 18 00 

RIGHT ASCENSION (J2000) 



Fig. 5. Same as Fig.|4] 



Mtot(< r) = 



gas' c 



where Z is the metallicity of the gas and fmetais,Q is the metal 
mass fraction of the Sun. The results are shown in Table 4. 



where x = r/rc, no is the central electron density, ji is the mean 
molecular weight in atomic mass (=0.6), G is the gravitational 
constant and nip is the proton mass. The best-fit parameters ob- 
tained from the spectral and spatial analysis are shown in Table 
4. We then evaluate Ra as the radius encompassing a fixed den- 
sity contrast with respect to the critical density p^. This is nec- 
essary to compare different clusters. Since Centaurus, is at very 
low redshift with the analyzed XMM-Newton observation we are 
looking at an overdensity of 4500. Thus we used such overden- 
sity to estimate the mass of the clusters. For the other four clus- 
ters we estimated also the parameter at an overdensity of 2500 
and at an overdensity corresponding to the area covered by the 
metal maps. We compute than the metal mass as : 



5. 1. Metal mass estimation using metallicity maps 

In the ionized intra-cluster plasma the ratio between the proton 
density n^ and the electron density n^ is approximately 0.82. 
Thus, the emission integral (EI) could be written as: 



J ripne 



EI = I nnnpdV ^ npngV ~ \.2npV. 



(5) 



For each spectrum the emission integral (EI) can be derived eas- 
ily using the normalization K of the thermal spectrum measured 
within XSPEC: 



^metals ^gas'^ f metal s,Q 



(4) 



EI = Kx\0'\Andl^g{\+zf]. 



(6) 
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Fig. 6. Significance maps of all the clusters. 
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Fig. 7. Plot of abundance against temperature for each bin. 
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Table 4. Best-fit results of the spectral and spatial analysis of the sample of galaxy clusters at an overdensity corresponding to 4500. 
We determined the total mass using the T derived fitting the spectra with a single temperature model. The last column represents 
the ratio between column 10 and 9 containing the metal mass. 



/?-model metal maps 



Cluster 




Z/Zo 






A 


^tot 


^^gas 


^^metals 


^^metals 


ratio 




keV 




kpc 






10^^ Mq 


10^^ M© 


10 M© 


10^^ M© 


C. 10 to C. 9 


Sersic 


2.38+0.01 


0.357+0.003 


37+1 


0.602+0.001 


4500 


2.59+0.18 


2.45+0.12 


1.75+0.09 


1.92+0.25 


1.34+0.23 


Centaurus 


3.74+0.01 


0.810+0.003 


17+1 


0.397+0.001 


4500 


2.91+0.14 


1.20+0.11 


1.94+0.19 


2.17+0.18 


1.12+0.19 


Hydra 


3.57+0.01 


0.357+0.005 


35+1 


0.561+0.001 


4500 


3.76+0.22 


3.36+0.20 


2.40+0.14 


3.32+0.59 


1.38+0.31 


A496 


3.72+0.01 


0.586+0.004 


27+1 


0.459+0.001 


4500 


3.76+0.19 


2.72+0.16 


3.19+0.18 


3.91+0.45 


1.23+0.20 


A2029 


7.01+0.02 


0.493+0.006 


41 + 1 


0.563+0.001 


4500 


12.26+0.81 


11.42+0.83 


11.79+1.01 


13.97+0.95 


1.18+0.16 


Sersic 


2.39+0.01 


0.309+0.003 


37+1 


0.602+0.001 


2500 


3.95+0.11 


3.95+0.21 


2.44+0.13 


3.03+0.21 


1.24+0.19 


Hydra 


3.46+0.01 


0.346+0.005 


35+1 


0.561+0.001 


2500 


5.47+0.29 


5.53+0.41 


3.83+0.43 


4.97+0.55 


1.30+0.26 


A496 


3.80+0.01 


0.506+0.004 


27+1 


0.459+0.001 


2500 


5.30+0.20 


4.79+0.31 


4.85+0.31 


6.38+0.70 


1.32+0.22 


A2029 


7.07+0.02 


0.468+0.006 


41 + 1 


0.563+0.001 


2500 


17.5+0.90 


21.8+1.30 


20.4+1.20 


23.7+2.80 


1.16+0.19 


A496 


4.43+0.01 


0.486+0.004 


27+1 


0.459+0.001 


1600 


7.78+0.36 


6.79+0.33 


6.60+0.38 


8.79+1.00 


1.33+0.21 


Hydra 


3.64+0.01 


0.299+0.005 


35+1 


0.561+0.001 


1400 


8.35+0.44 


10.29+0.56 


6.15+0.34 


7.63+0.99 


1.24+0.22 


Sersic 


2.31+0.01 


0.286+0.003 


37+1 


0.602+0.001 


1200 


5.05+0.34 


6.79+0.29 


3.88+0.29 


4.83+0.32 


1.24+0.16 


A2029 


7.81+0.02 


0.420+0.005 


41 + 1 


0.563+0.001 


1200 


34.0+1.80 


25.3+1.30 


21.3+1.04 


49.3+9.90 


2.31+0.55 



For each pixel the gas mass along the line of sight is determined 
using: 

M^as = ^ miUiW ^ (munu + rnHenHe)^ ~ 1.3m// Vv (7) 

where uu and uue are the proton and helium number den- 
sity respectively, mu is the proton mass, miie-^^p and V is 
the volume of the emitting region. This was determined (as- 
suming that the properties of the material in each region are 
constant and that there is no material projected onto them) as 
y ^ 2 V7?2 - X2 - Y^A, where A is the area of the region, R is 
the radius emcompassing the fixed density contrast, and X and 
Y are the projected distances in the east- west and north- south 
directions, respectively. We assumed the solar HejH fraction, 
nHelnn ~ 0.095 and we did not consider the mass contribution 
of ions heavier than He which are negligible compared to H and 
He. Then using the equation [4] we determined the metal mass 
along the line of sight for each pixel and summed them up for all 
the pixels. 

We found that the metal mass is higher than the metal mass 
obtained by assuming the metallicity mean of the investigated 
area (see Table 4). The explanation is that when we determine 
the metal mass using the mean metallicity we do not take into 
account properly that there are a lot of metals at positions of 
low density. Since, the metallicity is not constant throughout the 
galaxy cluster, as shown with the metallicity maps, using the 
maps we are estimating the metal mass along the line of sight 
better. The discrepancy changes from cluster to cluster, and con- 
sidering different radii. In the very inner part, at an overdensity 
of 4500 the metal mass can be underestimated up to 30%. At 
larger radii the discrepancies can be up to more than 2 times. 
From simulations we know that ram-pressure is more important 
than galactic winds in the center. The interstellar material in a 



galaxy feels the ram-pressure of the intracluster medium as it 
flows past. This ram-pressure is 

Pr^PlCMV^, (8) 

where p/cM is the ICM density and v is the relative velocity be- 
tween the galaxy and the ICM. 

From the virial theorem we know that the velocity of the galaxies 
is related to the total mass of the cluster. A higher mass of cluster 
will translate to a higher velocity of galaxy and due to the Eq. [8] 
to a higher metal mass as a consequence of ram-pressure strip- 
ping. iDe^GrandTetal] (12004 ) found that the iron mass associated 
with the abundance excess does not favor a scenario where the 
iron mass is accreted from the cooling flow and that the excess 
can be entirely produced by the brightest cluster galaxy (BCG) at 
the centre of cool core clusters. On the other hand in our sample 
of 5 clusters we found that in the very central part, at an over- 
density of 4500, the metal mass seems to be correlated with the 
total mass. The clusters with higher mass are able to strip more 
gas from the galaxies and to explain at least part of the central 
mass. 

6. Conclusions 

Based on XMM-Newton observations, we studied the spatial dis- 
tribution of metal abundances in a sample of 5 relaxed clusters. 
Below we summarize the main results. 

- Even for relaxed clusters the distribution of metals is clearly 
non- spherical. It looks very inhomogeneous with several 
maxima separated by low metallicity regions. 

- The deviation from the expected temperature-metallicity re- 
lation suggests that several processes are at work in galaxy 
clusters and that the simple picture of stripped gas does not 
hold. 



12 



L. Lovisari et al.: Inhomogeneous Metal Distribution in the Intra-Cluster Medium 



- The radial profiles is consistent with the current idea that the 
relative contribution of SN la to SNCC increases towards the 
cluster center and the 0/Fe is more sensitive to this ratio than 
Si/Fe. 

- Varying from cluster to cluster, the relative number of core- 
collapse supernovae necessary to reproduce the observed 
abundances ranges between 65-80% 

- Using a single extraction region to determine the metallicity 
gives a systematic underestimation of the metal mass: the 
metal masses are typically understimated by 10-30%. 
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